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ABSTRACT 

Context. The current generation of X-ray observatories like Chandra allows studies with very fine spatial details. It is now possible to 
resolve X~ray point sources projected into the cluster diffuse emission and exclude them from the analysis to estimate the "correct" 
X-ray observables. 

Aims. We wish to verify the incidence of point sources on the cluster thermal emission and to evaluate the impact of their non-thermal 
emission on the determination of cluster properties. 

Methods. To these ends we use a sample of 18 high-z (0.25<z<1.01) clusters from the Chandra archive and subtract the non-thermal 
emission of the point sources from the extended thermal emission due to the cluster itself. We perform a detailed analysis of the 
cluster properties and compare the changes observed in the X~ray observables, like temperature and luminosity or their inter-relation, 
when one keeps the point sources in the analysis. 

Results. The point sources projected into the cluster extended emission affect the estimates of cluster temperature or luminosity 
considerably (up to 13% and 17% respectively). These percentages become even larger for clusters with z>0.7 where temperature and 
luminosity increase up to 24% and 22%, respectively. 

Conclusions. The conclusions are that point sources should be removed to correctly estimate the cluster properties. However the 
inclusion of the point sources does not impact significantly the slope and normalization of the L^oj-T relationship since for each 
cluster the correction to be applied to T and Lj,„; produces a moderate shift in the L(,o;-T plane almost parallel to the best-fit of the 
"correct" h/^oi-T relation. 
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1. Introduction 

O 



X-ray studies of clusters of galaxies performed with X-ray teles copes like ASCA (Tanaka et anil994l) . Beppo-SAX dParmar et al.l 
119971: iManzo et al.ll 19971: iBoella et alJl 19971: iFrontera et anil997[) or ROSAT ( lTrumpeiill984h have attributed the total X-ray emis- 



sion from clusters of galaxies to thermal bremsstrahlung emission from the thin hot gas that fills the regions between the cluster 
galaxies. This is due to the optimal consistence between the observed X-ray spectra and the expected thermal emission from the 
5^ highly ionized hydrogen and helium in the intra-cluster medium at temperatures T ~ 10^ - 1.5 x 10^ K. The detection of point 
sources in the cluster fields has been hindered b y the traditionally poor angu lar resolution of su ch X-ray telescopes . The current 
gene ration of X-ray observatories like Chandra dvan Spevbroeck et alj|199'7l) or XMM-Newton dStriider et al.ll200TI : [Turner et al.l 
12001 ) has revolutionized X-ray astronomy, enabling studies with very fine spatial details. In addition the improvement of detection 
techniques (e.g. the multiscale wavelet detection bv iFreenian et all 1200 21) can now reliably separate small-scale source emission 
from surrounding larger-scale diffuse cluster emission. Hence, it is now possible to subtract the non-thermal emission of the point 
sources embedded in the cluster emission from the extended thermal emission due to the cluster itself. In recent years a number of 
detailed XMM-Newton and Chandra studies have allowed investigators to compute the "correct X-ray observables" and thus the 
"correct scaling relations" between them excluding the point sources from the analysis of the cluster emission. 

In this paper we present a detailed analysis of the archival Chandra data for a sample of high-z clusters. The aim is to verify the 
incidence of point sources on the thermal emission of the cluster and to evaluate their contribution to the cluster total emission. In 
particular we will examine the point source effects on the cluster temperature and luminosity and on t he L^-T relation . The sample 
analyzed consists of eighteen clusters with redshift in the range 0.25 < z < 1.01, that was also used bylBranchesi et al.l (12007) (from 
now on BR07) to check for any overdensity of point sources in the inner region of clusters of galaxies. In a companion paper (Paper 
II: Branchesi et al. 2007) we will use these same clusters in combination with clusters taken from the literature to revisit the "correct 
Lx-T" relation and check its evolution with redshift. 

All uncertainties in this work are at the 1 cr confidence level, unless otherwise noted. We use a ACDM cosmology with Ho = 
70 kms ' Mpc ' and Q,n = 1 - = 0.3. 
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Table 1. Cluster sample parameters and details of Chandra observations 



Cluster name 




z 


RA 


DEC 


Obs.ID 


ACIS 


Mode 


Exp. 


Nh 








hh 


mm ss 






II 








ks 


102« 


cm - 


(1) 




(2) 




(3) 




(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


Abell2125 





246 


15 


41 


12 


+66 


16 


01 


2207 


I 


VF 


79.7 


2 


77 


ZWCL 1454.8+2233 





258 


14 


57 


15 


+22 


20 


33 


4192 


I 


VF 


91.4 


3 


22 


MS 1008.1-1224 





302 


10 


10 


32 


-12 


39 


23 


926 


I 


VF 


44.2 


6 


74 


ZWCL 0024.0+ 1652 





394 


00 


26 


35 


+ 17 


09 


39 


929 




VF 


36.7 


4 


19 


MS 1621.5+2640 





426 


16 


23 


36 


+26 


34 


21 


546 


I 


F 


30.0 


3 


59 


RXJ 1701.3+6414 





453 


17 


01 


24 


+64 


14 


10 


547 


I 


VF 


49.5 


2 


59 


CL 1641+4001 





464 


16 


41 


53 


+40 


01 


46 


3575 


I 


VF 


44.0 


1 


02 


V 1524.6+0957 





516 


15 


24 


40 


+09 


57 


48 


1664 


I 


VF 


49.9 


2 


92 


MS 045 1.6-0305 





539 


04 


54 


12 


-03 


00 


53 


902 




F 


41.5 


5 


18 


V 1121+2327 





562 


11 


20 


57 


+23 


26 


27 


1660 


I 


VF 


66.9 


1 


30 


MS 2053.7-0449 





583 


20 


56 


21 


-04 


37 


51 


1667 




VF 


43.5 


4 


96 


V 1221+4918 





700 


12 


21 


26 


+49 


18 


30 


1662 




VF 


79.4 


1 


44 


MS 1137.5+6625 





782 


11 


40 


22 


+66 


08 


18 


536 




VF 


117.5 


1 


18 


RDCSJ 1317+2911 





805 


13 


17 


21 


+29 


11 


19 


2228 




VF 


111.3 


1 


04 


RDCSJ 1350+6007 





805 


13 


50 


48 


+60 


06 


54 


2229 




VF 


58.3 


1 


76 


RXJ 1716.4+6708 





813 


17 


16 


49 


+67 


08 


26 


548 




F 


51.5 


3 


71 


MS 1054.4-0321 





830 


10 


56 


59 


-03 


37 


37 


512 


S 


F 


67.5 


3 


67 


WARPJ 1415.1+3612 


1 


013 


14 


15 


11 


+36 


12 


00 


4163 


I 


VF 


89.2 


1 


10 



- Column 1 : Cluster name 

- Column 2: Spectroscopic redshift tabulated in the literature 

- Column 3-4: Right ascension and declination (Equatorial J2000, HH MM SS.S, +DD MM SS.S) of the centroid of the Chandra photon 
distribution in the 0.5-5 keV energy band assumed as the cluster center 

- Column 5: Identification number of the observation 

- Column 6: Detector where the aimpoint lies (I, for ACIS-I or S, for AC/5-5) 

- Column 7: Observation mode (F for FAINT or VF for VFAINT) 

- Column 8: Exposure time in ks corresponding to the nominal exposure filtered to exclude time periods of high background 

- Column 9: Column density of Galactic hydrogen in units of lO"" cm"^, obtained from the Chandra X-ray Center (CXC) Proposal Planning 
Tool Colden (Galactic Neutral Hydrogen Density Calculator): NRAO-compilation by Dickey & Lockman ( 1990) 



2. Chandra data reduction and cluster analysis 

We extracted from the Chandra archive the X-ray data of 18 galaxy clusters with redshift in the range 0.25< z <1.01 and exposure 
time greater than 30 ks. Data reduction was performed using version 3.2.1 of the CIAO software (Chandra Interactive Analysis of 
Observations; see web page http://cxc.harvard.edu/ciao/index.html) and version 3.0.3 of the CALDB (Calibration Database). The 
sample is the same used by BR07 where the selection criteria and the data reduction are described in more detail. For clarification 
purposes we report here Table[T]of BR07 that lists the sample parameters and the details of the Chandra observations. Clusters are 
arranged in increasing redshift order. 

The clusters of galaxies analyzed here appear as extended sources in the Chandra images. Differently from lower resolution instru- 
ments, Chandra is able to reveal point sources overimposed but not necessarily associated with the pointed cluster. Although such 
objects may be scientifically interesting in their own right (see discussion in BR07) investigators can now filter out any point source 
before fitting spectral models to the extended emission. For a correct analysis of the cluster properties the point source (non-thermal) 
emission should be subtracted from the cluster extended emission. Since one of the goals of our analysis is to estimate the impact 
of these point sources on the cluster properties, two approaches have been adopted: 

- Point sources are identified using the CIAO Detect package WAVDETECT, which has the ability to work in complex fields with 
both point and extended sources (see BR07 for details on the procedure followed). The point sources are later removed from the 
event file using the dmcopy command. 

- Point sources are not removed from the event file in an attempt to evaluate their efl'ects on cluster parameters. 

2.1. Background subtraction 

The issue of background subtraction was carefully considered. For both spatial and spectral analyses, it is necessary to associate 
a background to the source events. For most of the clusters the background was measured locally, within the same target field, 
in a region free of point sources and in the vicinity of the cluster but not contaminated by the cluster emission. We checked that 
variations of the background intensity across the chip do not affect the background subtraction by comparing the count rates in the 
cluster and in the background regions at energies larger than 8 ke V, where the signal from the cluster should be nil. 
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We also used the "blank-sky" data sets, following the steps and recommendations of the Markevitch's COOKBOOKQ. The "blank- 
sky" background files were first reprocessed and reprojected to match the corresponding cluster observation gain and position. 
The background files were then normalized to the shorter exposure time of the cluster observations. The background events were 
extracted from the same region of the chip as the cluster to model spatial variations in the background. Obtaining background spectra 
from blank-sky data sets has the advantage to use the same region as the source, thus eliminating potential systematic errors caused 
by spatial variations of both the energy response and the effective area across the chip. The method remains, however, vulnerable to 
temporal variations in the spectrum of the particle background and also cannot easily account for the strong directional variation of 
the Galactic soft X-ray emission. 

We performed consistency checks for several clusters of the sample. The properties obtained using background files derived from 
different regions of the target field and from the "blank-sky" data sets are generally independent of the background used. For the 
two most nearby clusters of the sample, ZWCL 1454.8-H2233 and Abell2125, which cover a very extended region, we considered 
more appropriate to extract the background spectrum file from the "blank-sky" data-sets. For all the other clusters the background 
spectrum file was extracted from regions within the same target field. 



2.2. Spatial analysis 

2.2.1. Cluster surface brightness 

The spatial analysis of the cluster X-ray emission was performed within images (with point sources excluded) in the energy range 
0.5-5.0 keV. Under the assumption that th e cluster gas is a symmetric isothermal sphere with a density profile described by a 
jS- model dCavaliere & Fusco-Femianoll 19761) : 



UgasCr) - no,gas 



-3/3/2 



(1) 



the radial surface brightness was modelled accordingly as: 



S(r) = So 



2Y 



8+1/2 



(2) 



where Rc is the core radius and /3 describes the slope of the density profile at large radii. 

For each cluster the yS-model fit to the surface brightness was performed using the exposure-corrected image with a constant back- 
ground included in the fit. The best-fit parameters for Rc and /3 are listed in Table |2] 

Even if the simple jS-model is known to poorly describe the radial profile of Chandra highly resolved clusters (see lEttori et al.l2004h . 
this model provides a good description of the cluster for our pourpose, that is to extrapolate luminosities to larger radii (see Sect. [3] 
and Paper II of Branchesi et al. 2007). Exceptions are MS 1054.4-0321 and ZWCL 1454.8+2233. The former cluster presents a 
clear western substructure. The fit to the whole cluster provides unreasonably large values for /3 and Rc (as previously noted by 
I.Teltemaetal.l '200n. The fit improves once we mask the western substructure from the main body of the cluster. The values of Rc 
and yS in Table|2]for MS 1054.4-0321 were obtained in this way. For ZWCL 1454.8+2233 a single y8-model is not an acceptable 
representation of the radial profile (see Notes on individual clusters in Sect. |4|i. Therefore no values of Rc and /3 are indicated in 
Table m 



2.2.2. Definition of cluster extent 

The spectrum of each cluster was extracted from a circular region centered on the X-ray centroid, out to a maximum radius which is 
hereafter referred to as spectral radius, Rspec- This radius was chosen individually cluster by cluster so as to optimize the signal-to- 
noise ratio in order to determine the X-ray temperature of the cluster with maximum count statistics. To define this circular region 
around the centroid of the photon distribution we followed the procedure suggested bv lTozzi et al.l ( l2003h . 

A second extraction radius was adopted in order to consider the cluster extent where diffuse emission is still detectable. This radius, 
named Rext, is the radius where the cluster radial surface brightness merges into the background, and beyond which no further 
significant cluster emission is detected. Both Rspec and Rext are listed in Table|2] The fraction of the net counts included in the Rspec 
extraction region is always between 0.80 and 0.95 of the net counts included in the Rext extraction region once point sources have 
been excised from the cluster emission. 

For each cluster, the events included in the extraction region were used to produce a spectrum file. This procedure was repeated 
considering the two above defined radii (Rspec and Rext) and approaches, that is filtering out or not the events associated to the point 
sources from the event file. We ended up with four spectrum files per cluster. 



' See http://cxc.harvard.edii/contrib/inaxim/acisbg/data/ README\ ?(nd 
\http://cxc.harvard.edu/contrib/maxim/acisbg/COOKBOOK\ 
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2.3. Spectral analysis 

Spectra were extracted from within both the defined detection radii, Rspec and Rext, considering the two situations described above: 
cluster emission alone (from now on "cluster") and cluster plus point source emission (from now on "cluster+ps"). The Auxiliary 
Response File (ARF) and the Redistribution Matrix File (RMF) were computed from the same region where the spectra are extracted 
and were weighted by the detected counts in the restricted energy range (0.5-2.0 keV) where both the response and the thermal 
model do not vary much. When possible the weighted RMFs have been created using the new CIAO tool mkacismij^ otherwise 
the previous tool mkrmf'WcLS used. For each cluster spectrum file (with or without point sources) the respective ARFs and RMFs 
were generated in order to account for the presence or absence of point sources. For each cluster spectrum file there are then three 
associated files, namely the background spectrum file and the two response matrices. 

The spectra are then analyzed with the XSPEC package and fitted over the energy range 0.8-7.0 keV. The photons with energy below 
0.8 keV were excluded to avoid systematic biases in the temperature determination due to uncertainties in the ACIS calibration at low 
energies. Ignoring energies above 7 keV has little effect on Chandra data due to the low effective area above that energy and the rapid 
S/N decrease of the thermal spectra . The spectra were fitted with an absorbed single-temperature thermal model called wabs(mekal) 
(lKaastral l l992t lOed ahl et al.l ll995h. The absorbing hydrogen column was frozen at the Galactic value (as determined from radio 
HI maps, Dickey & Lockma^ 19901) in correspondence of the X-ray peak. The gas temperature, T, and the normalization, K, of 



the thermal component are the only free parameters. The best-fit temperatures were determined freezing the redshift at the values 
measured by spectroscopic observations (available in the literature) and fixing the metallicity (Z) at 0.3 Z©. In the spectral fitting, we 
take into account the increased effective area at energies larger than 2 keV (due to a thin hydrocarbon layer, see Marshall et al. 2003|) 
including in the fitting model the "positive absorption edge" (XSPEC model edge) described by Eq. 1 in iVikhlinin et al.. (2005) . 

The spectral fit has been performed using both the Cash statistics ( ICashll9'79b and the;^'^ statistics (adopting a standard binning with 
a minimum of 20 phot ons per energy chann el in the source plus background spectrum). The Cash statistics seems to be preferable 
for low-S/N spectra (Nouse k & Shuel[l989 *) when the number of counts available per bin is low. However the agreement between 
the models obtained with the two statistics (see Appendix lA.il ) makes us confident that both statistics are a good choice. Throughout 
this paper we will use the statistics for the only reason that, unlike the Cash statistics, the statistics gives a measure of the 
absolute goodness of the fit. We verified that all our spectra are well-fitted by single-temperature models: the best-fit models have 
a reduced ~ 1 and a null-hypothesis probability above 15% (except for ZW CL 1454.8+2233, see Notes on individual clusters in 
Sect.©. 

Once the models with the best-fit parameters were determined XSPEC has been used to calculate the cluster flux over different energy 
bands and the cluster bolometric luminosity. The errors on the cluster parameters are obtained in XSPEC from the distribution of 
the values around the best-fit value of the spectral analysis. The quoted luminosities were corrected for the effect of absorption by 
the Galactic HI column density at low energies. Hereafter the unabsorbed bolometric luminosity is called bolometric luminosity or 
Lboi- 

Since the quality of the fits is better inside the region that maximizes the signal to noise ratio, the temperatures estimated within 
Rspec, i-e- Tr^p.,, are considered more representative of the actual average temperature of the gas. However the two temperatures are 
consistent within the errors with no significant systematic differences between them. For the luminosities we used instead the Rext 
radius in order to take into account the faint brightness tails at the cluster periphery as well as the point sources in those regions. 
The choice of the different radii (Rspec for the determination of temperature and Rext for the determination of luminosity) is justified 
in Appendix lA.2l where a more detailed description of the analysis of the best-fit cluster parameters is given. 

The cluster parameters are listed in Table|2] The columns contain the following information: 

- Column 1 : The first line gives the cluster name, the second and third lines indicate how parameters are derived, if for the cluster 
alone or for the cluster pl us point sources . An as terisk close to the cluster name indicates that the cluster has been classified as 
a possible cooling core bv lvikhlininetan (120021) 

- Column 2-3: Core radius in kpc and /3 

- Column 4: Radius (spec) which maximizes the S/N, in arcsec and kpc 

- Column 5: Radius (ext) where the cluster X-ray radial profile becomes flat, in arcsec and kpc 

- Column 6: Number of point sources detected (and removed) within Rspec and Rext, respectively. These point sources were 
detected either in the soft or in the hard energy band (see BR07 for details). The column lists also faint sources which become 
significant (S/N > 3) in the full energy band (0.5-7.0 keV). 

- Column 7: Temperature estimated within Rspec in keV 

- Column 8: Temperature estimated within Rext in keV 

- Column 9-10: Observed soft and hard flux estimated within Rext in units of 10 erg cm"^ s ' 

- Column 11: Unabsorbed X-ray bolometric luminosity estimated within Rext in units of lO'*'* erg s ' 



^ mkacisrmfhdis, been calibrated for the ACIS-I array plus ACIS-S1,S2, and S3. The tools creates response files intended for use with -120° C 
data that has the time-dependent gain adjustment and CTI correction (if available) applied. There is no CTI correction for the back-illuminated 
ACIS chips, SI and S3. 
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5 



Cluster 
(1) 


Rc 
(2) 


P 

(3) 


Rspec 

(4) 


R-ext 

(" ) (kpc) 
(5) 


N 

^ ^ sources 
(Rspec) (Rext) 

(6) 


Tr 

*^spcc 

(keV) 
(7) 


TRex, 

(keV) 
(8) 


So.5-2.0 
('10-13 c?s~) 

(9) 


S2.0-10.0 

nO"13 Cffs") 

(10) 


Lfoo/ 

(11) 


Abell2125 

CI 

cl-t-ps 


182±12 


0.54+0.02 


153 591 


241 


931 


13 


26 


'1 4.+0.2 T c+0.2 
-'•^-0.2 -^-^-0.2 
^ c+O.T a 7+0.3 
-'-^-0.1 -'•'-0.2 


4 06+0-14 
4 70+0.14 
^■'^-0.14 


4 1Q+0-21 

5 

-'•^^-0.23 


9 9n+0.06 

2 64+0-0^ 

^•"^-0.07 


ZWCL 1454.8+2233 
cl 






128 512 


200 


800 


4 


7 


4 4+0.1 


4.5+01 
4 5+^-' 


32 05+0-18 


on 71 +0.44 


20.83!0.;o 


cl+ps 
















4 4+0.1 


^7 Q4+0.I7 


40 08+0-42 


21 02+0- '0 


MS 1008.1-1224 
cl 


165±9 


0.64+0.02 


128 572 


172 


771 


6 


10 


6 0+"'* 
-Q-3 


6 1+0-4 
6 3+0-4 


8.76+0-iJ 
8 97+0-20 


15 67+0-«' 

16 32+0-''3 


10.2810-29 


cl+ps 
















6 2*^-^ 

"•^-0.4 


10.65+0-30 


ZWCL 0024.0+1652* 
cl 

cl+ps 


128±10 


0.67+0.03 


69 367 


118 


628 


2 


7 


4 4+0.5 
4 8+° *' 

^•°-0.6 


4 4+0.7 

5 0+^-^ 

-'•^-0.7 


9 9C+0.14 
^■^-'-0.14 

2 50+0-12 

^■-'"-0.15 


2 67+0-25 

^•"'-0.28 

T TC+0.33 

•'•-'-'-0.40 


3 90+0-22 

-'•^"-0.22 

4 55+0.26 
^•-'-'-0.29 


MS 1621.5+2640 
cl 

cl+ps 


227±17 


0.65+0.03 


118 659 


148 


823 


3 


5 




7 5+1-3 

^•0!|:^ 


4 T 1+0.21 


7 -70 +0.41 
'•'^-0.59 


1 1 07+0-54 
' '-"'-0.56 

1 ') 69+0-™ 

'-•"^-0.82 


RXJ 1701.3+6414* 

Cl 

cl+ps 


15±2 


0.41+0.01 


79 455 


108 


626 


2 


3 


4 Q+0.5 c c+0.5 
•°-n.4 -'"'-O.fi 


2.61+0-14 
2 62+^-i^ 

^•"^-0 15 


o 9A+0.24 

3 55+0-^* 


97+0.28 

6 54+0-28 

"••^ -0 37 


CL 1641+4001 
cl 

cl+ps 


151±18 


0.77+0.06 


54 317 


89 


519 


3 


6 


5 1+0.8 5 5+1.1 

4 Q+0.6 5 7+0.9 
^•^-0.6 •'■'-0.7 


1 05+0-10 
1 57+0-12 

^•-"-0.13 


1 36+0-14 

2 09+0-15 

^■"^-0.19 


2 65+0-22 
4 02+0-22 

^•"^-0.26 


V 1524.6+0957 
cl 


302±27 


0.80+0.05 


79 488 


148 


916 


3 


10 




2 00+0-18 


2.7210-22 


6.85+0-44 


cl+ps 
















-'•^-0.6 


6 4+'-' 


9 OQ+0.I6 
^•-'^-0.17 


3 57+0-32 


8 58+0-50 

°-'°-0.64 


MS 045 1.6-0305 

Cl 

cl+ps 


270±8 


0.90+0.02 


89 562 


148 


937 


3 


6 


9-4-05 
9 3+o:5 


9.6+1-0 
9 8+*^-^ 

^•°-0.8 


q 81+0-20 

-Q?Q 
9 96+0-19 

='•='"-0.24 


20.84:[-i2 
21 41+0-'9 

■ -1-44 


Sfl 4^+2-56 
JU.M-j_2 

51 77+2:44 
■"-2.44 


V 1121+2327 
cl 

cl+ps 


437±58 


1.19+0.18 


67 434 


128 


829 


4 


8 


4 c+0.5 c c + 1.1 
^•-'-0.4 -'"'-0.9 
4 "n+O-S c q+1.2 


1 58+o-;3 


1 79+0.19 
^•'^-0.23 

2 08+0-12 

■""-0.26 


5 44+0.37 
-'•^^-0.43 

6 38+0-35 

"•-^"-0.46 


MS 2053.7-0449* 
cl 


115+12 


0.64+0.03 


57 373 


118 


779 


1 


3 


4 3+0-5 
-Q-4 


5 1+1-4 

-'•^-1.0 


1 26+0-12 


1 51+0-19 
-Q-33 
1 56+0-20 


5 59+0.46 

-Q-§i 


cl+ps 
















4 5+0.6 
^•-'-0.4 


5 2+1-4 


1 28+0-;^ 


5 81+0-43 

-0.60 


V 1221+4918 

Cl 

cl+ps 


272+20 


0.76+0.04 


79 562 


143 


1020 


3 


8 


7.0!0.8 6.4:0.9 

7 9+0.8 7 o+l.l 
'• -0.3 '-^-0.8 


1 qn+0.10 

9 1Q+0.11 
^•^=^-0.12 
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Fig. 1. Comparison between the estimates of the gas temperature by VI02 {left panel) and our estimate for the 1 3 clusters in common. 
To the right panel the same comparison is shown for the 12 clusters in common with ET04. The circled points indicate clusters 
assumed by VI02 to be cooling core clusters. Like us, ET04 did not excise any excess surface brightness central region due to the 
presence of "cooling flows". In both panels the dashed line is equality between the two works. 

3. Comparison with other authors 

Two recent works, one by'Vikhlinin etalj ( l2002 h (from now on VI02) and one bv lEttori et al.l (|2004') (from now on ET04), have 13 
and 12 clusters respectively in common with us. It is thus very instructive to compare our estimates of the gas temperatures and 
bolometric luminosities with their results in order to check for the presence of any systematic bias due to the different approaches 
adopted. Of particular importance is the radius used to measure the properties of the clusters. 

3.1. Temperatures 

A comparison of our temperatures with those of VI02 is shown in Fig. [T| (panel to the left). Differently from us, VI02 excluded 
the central 100 Ii^q kpc region in the cooling core clusters. The agreement is very good: a mean ratio of 1.02+0.04 between VI02 
temperatures and ours has been found. In three out of four clusters assumed by VI02 to be possible cooling core (indicated by an 
asterisk in Table |2]) they report higher temperatures. The largest difference, significant at 2.9cr, is found for RXJ 1701+6414. For 
the fourth cluster, RDCSJ 1317+291 1, we find instead a higher temperature, kT - 3.7^^^ keV against the value found by VI02 of 
kT = 2.2;^°;^ keV. The difference is significant at the 2.2cr level. For MS 0451.6-0305 the temperature difference, significant at 2cr 
confidence level, can be explained in terms of the new calibrations that we have applied. 

The comparison with ET04 (their values vs. ours) shows that their temperatures are on average a factor of 1 .06+0.03 higher than ours 
(see panel to the right in Fig.lTJ- The largest absolute values of the differences (significant at > 2cr) are found for MS 0451.6-0305 
(3.6cr), MS 1054-0321 (2.5cr) and MS 2053.7-0449 (2.2cr). All these discrepancies are discussed in the notes on individual clusters 
(Sect. Si. 

3.2. Luminosities 

It has to be noted that both VI02 and ET04 derived cluster luminosities within areas often quite different from the areas we used. 
Many authors discuss the importance of the choice of the radius within which cluster properties are measured, especially when 
comparing integrated cluster properties with theoretical predictions or simulations. A common choice is to use as radius a fixed 
linear size, which has the obvious benefits in terms of simplicity. VI02 used a fixed radius of 2 h^^^ Mpc, corresponding to 1 .4 /i^g 
Mpc in our cosmology. Since VI02 excluded the central 100 /i^g kpc regions in the cooling core clusters, they accounted for the 
missed flux by multiplying by a factor 1.06 typical of a y6-model cluster 

Other authors define a "more physical" radius which requires the knowledge of the cluster mass profile so that the mean enclosed 
density is a fixed factor above the critical density of the Universe. This approach is followed by ET04 who define R500, which 
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Fig. 2. Comparison between the luminosities estimated by VI02 and by us for the 13 clusters in common (left panel) and by ET04 
and by us for the 12 clusters in common (right panel). The circled points indicate the clusters assumed by VI02 to be cooling core 
clusters. In both panels the dashed line is equality between the two works. 



corresponds to an overdensity A, - 500 x A,,(z)/(18;r^) with respect to the critical density of the Universe at redshift z. Their 
luminosities were computed by extrapolating to R500 the values measured within Rspec by means of an isothermal j6-profile. Like 
us, ET04 did not excise any excess surface brightness central region due to the presence of "cooling flows". Differently from these 
authors we used an "observed" radius which indicates the region where the emission is detected (Rext, see Sect. 12.2.21 ). The VI02 
radii differ from our radii by a factor ranging from 1.04 up to 2.7. Also the ET04 radii are larger than our radii by a factor ranging 
from 1.03 up to 1.55, except for RDCSJ 1350+6007 that has a ET04 radius about 30% less than ours. 

The comparison between the VI02 luminosities and ours gives a mean ratio of 1.10+0.04 (see left panel of Fig.|2]l. The fact that the 
VI02 luminosities are on average 10% higher than our estimates may be explained by the larger regions (radius of 2 /j^y Mpc) used 
by them. The mean ratio between the ET04 luminosities and ours is 0.96+0.03 (see right panel of Fig.|2|. 

To check if any systematic bias is introduced by our choice of the radius, we extrapolate our luminosities to the radii adopted by the 
other authors. The correction to be applied to our data was obtained assuming that the gas density profile is described by the y6-model 
(see Eq.[T]) defined by our estimates for the core radius, Rc, and for /3 (see Sect. 12. 2. Tl ). It has to be noted that these corrections do not 
depend on the ratio between the different integration radii only but depend strongly on and Rc- This can be clearly seen in Fig. [3] 
which shows the luminosity computed within a given radius R as a function of R/R^, for different values of /3. 

The comparison between VI02 luminosities and our luminosities extrapolated to 1 .4 h^y Mpc gives a very good agreement (see 
left panel of Fig.lUi. Excluding the most discrepant cluster in the plot, RDCSJ 1317+2911 (which has although very large errorbar 
values) the mean ratio between the VI02 luminosities and ours is of 0.99+0.03. Note that the corrections to be applied to our 
luminosities for such an extrapolation are not negligible. They range from ~ 0.9 to ~ 1.25 with a median value of 1.10 (mean 1.13). 
Similar results were obtained using the VI02 estimates of and /3, instead of ours, to extrapolate to the VI02 radius. The clusters 
with the largest absolute values (> 2cr) of the difference between our and VI02 luminosities are indicated in the figure: the possible 
cooling core cluster RXJ 1701.3+6414 (4.0o-), MS 1 137.5+6625 (3.2cr) and MS 0451.6-0305 (2.3cr). 

In the comparison with ET04 (see right panel in Fig.lUi we extrapolated our luminosities to R500. This radius was calculated assuming 
our estimates for the cluster temperature, for the core radius and /3 (see Paper II of Branchesi et al. 2007 for more details). If one 
compares the right panel of Fig. |4] with the right panel of Fig. |2] it can be noted that the extrapolation to R500 does not improve the 
comparison with ET04. In fact the mean ratio between the ET04 luminosities and ours changes from 0.96+0.03 to 0.93+0.02. A 
similar result is obtained if we extrapolate our luminosities using their estimates for R500, Rc and /3. In the right panel in Fig.|4]cluster 
luminosities which deviate more than > 2cr from the ET04 estimates are indicated with the cluster name. They are MS 1621.5+2640 
(4.0cr), MS 1054.4-0321 (3.0o-), MS 1 137.5+6625 (2.7cr). 

In summary the comparison with VI02 suggests that the correction for the different radii is important if we want to be consistent 
with their luminosities. Their radius is in fact quite larger than ours, but after corrections are implemented the agreement is very 
good. The comparison with ET04 instead suggests that the correction for the different radii are smaller than the measurement 
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Fig. 3. The luminosity computed within a radius R as a function of the ratio between the radius R and the core radius Rc for different 
j6 values. 

uncertainties and hence could be neglected. The agreement between our and ET04 measurements is good, although there is a very 
weak systematic offset in the sense that ET04 luminosities tend to be lower than our estimates. 

In paper II of Branchesi et al. 2007 we describe in detail why it is important to extrapolate all luminosities to an homogeneous radius 
(e.g., R500) when the observed Li,o/-T relations are compared with the self-similarity evolution predictions. 

4. Notes on individual clusters 

4.7. MS 0451. 6-0305 

MS 0451.6-0305 is the most luminous cluster in the EMSS sample (iGioia et alJll990h . We found a best-fit temperature within Rspec 
of kT = 9.4+°^ keV, which disagrees with the estimates of kT = 8.1+°-^ keV and of kT = 8.0!°-^ keV found by VI02 and ET04, 
respectively (see Fig.lTJ. 

In a more recent article lOonahue et al.l (l2003h analyzed the same Chandra data. They discuss in detail the results of applying a soft- 
energy, time-dependent correction to the ACIS-S, which however they consider uncertain. Thus without applying the correction the 
authors find that MS 045 1 .6-0305 is consistent with an isothermal cluster with kT ranging from 10 keV to 10.6 keV (± 1 .6 keV at the 
90% confidence level), and with intracluster Fe abundance range between 0.32 and 0.40 (+ 0. 13 solar at the 90% confidence level). 
Including the correction in their analysis, they may explain the discrepancy between their best-fit temperature and the temperature 
obtained by VI02. However they find that to be acceptable the fit requires a second component that could be either a cooler thermal 
component or a steep power-law component. 

Our data have been analyzed applying the time-dependent coiTection as suggested by the C XC (Chandra X-ray C enter). Our 
correction is more accurate since the CALDB used by us is more recent than the one used by Donah ue et alj (|2003). The more 
recent calibration adopted can explain the temperature discrepancies with resp ect to VI02 and ET04 (see F ig.[T|). Our estimate of 
kT = 9.4+IJ^ keV is consistent with the ASCA measure of kT = 10.2+[| keV dMushotzkv & ScharCl997l) . We detected six very 
faint sources (So.5-1 okev < lO"'"* erg cm"^ s"') within the more extended aperture radius Rext = 100" (see BR07). The cluster 
temperature does not change much if one includes the point sources (kT - 9.6^q ° keV vs kT = 9.8^p g keV). 

4.2. MS 1054.4-0321 

MS 1054-0321 is th e highest redshift ( z=0.83) cluster in the EMSS and shows a significant amount of substructure with the 
Chandra resolution dJeltema et al.ll200Tl) . As mentioned in Sect. [3] the temperature and luminosity estimates obtained by ET04 
disagree with ours (see right panels of Fig. [T] and Fig. H] respectively). ET04 estimated both the temperature and the best-fit sur- 
face brightness profile from the main body of the cluster after masking with a 36" radius circle the cooler region at RA, Dec 
(2000)= 10''56™55^7, -03°37'37". Since we did not exclude this region a lower temperature and a higher luminosity are obviously 
estimated. In addition, this cluster is an example of how t he X-ray ternperatu res based on Chandra data change as new calibrations 
become available. An analysis of XMM-Newton data by iGioia et al.l d2004h results in a temperature kT - 7.2^^^ keV, which is 
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Fig. 4. {Left panel) Comparison between the luminosities estimated by VI02 and our luminosities extrapolated to 1 .4 h^^ Mpc for 
the 13 clusters in common. The circled points indicate the clusters assumed by VI02 to be cooling core clusters. {Right panel) 
Comparison between the luminosities estimated by ET04 and our luminosities extrapolated to Rsoo for the clusters in common. In 
both panels the uncertainties on the extrapolated luminosities are combined errors which take into account both the errors obtained 
in the spectral analysis by XSPEC and the uncertainties on the y6-model for the gas density profile. In both panels the dashed line is 
equality between the two works. 

much lower than the temperature previously reported from A5CA data, kT - 12.3^^ 2 keV ("Don ahue et al.|[T998l) . and also somewhat 
lower than the first C handra temperature, kT - 10.4^[ ^ keV, determined by Jeltema et al.. (,2(joTl) . The temperature measurement 
of MS 1054-0321 by I Jeltema et al.l (1200 ll) probably suffered from the absence of a low-energy correction, called ACISABS, which 
was not available at the time of their analysis. VI02 used the same Chandra observations and derived a lower value for the tempera- 
ture, kT - 7.8 + 0.6 keV, in agreement with us and with the determination by Tozzi et al. (2003) of kT = 8.0 + 0.5 keV. All quoted 
uncertainties are at 90% confidence level except for the last one which is at 68%. 

As described in section [23] we obtain an estimate for the temperature of ^T - 8.3^" ^ keV within Rspec = 80". We applied 
the new available procedure to correct for the quantum efficiency. As observed in Jee et al. (2005) this new procedure (certainly 
more accurate) tends to yield a higher temperature than the ACISABS prescription. iJee et al.l (l2005b find a temperature of kT = 
8.9;^Q g keV, within 90". Using the more extended region {Rexr- 130") we estimate a temperature of kT = 7.8;^j q keV. 

4.3. MS 1137.5+6625 

MS 1137. 5-1-6625 is the second most distant cluster in the EMSS sample. Our best-fit temperature of kT - 6.2^^ g keV is consistent 
with the estimates obtained by VI02 (kT = 6.3+°^ keV) and by iBorpani et al.1 ( l200lb (kT = 5.7+jJ^ keV). The three temperatures 
above are computed within an ap erture radius o f abou t 60". Our temperat ure i s also consistent w ith the one determined from 
ASCA data, kT = 5.7+J}^ keV, by iDonahue et all (Il999l) . lEttori et all (|2004|) and lTozzi et all (1200 ll) found kT = 6.9!°-^ keV and 
kT = 7.0^^ 5 keV, respectively, within a smaller region of about 50" radius. 

4.4. RDCSJ 131 7+2911 

Despite this cluster is classified by VI02 as a possible cooling core system, our temperature estimate (kT = 3.7^^^ keV) is higher 
than the value found by VI02 (kT - 2.2;^^^ keV) even though we did not exclude the cooling flow region in our analysis. This 
discrepancy, visible in Fig.[T](left panel), might be explained considering that RDCSJ 131 7-H291 1 has a low signal to noise ratio. 
A discrepancy in the same direction has been found by ET04 (kT - 4.1^Q g keV) and by iTozzi et al.l (l2003b (kT - 4.0^^ ^ keV). 
iTozzi et al.1 (l2003l) argue that such a difference can be ascribed to differences in the procedure used to remove faint point sources 
within the extraction region, which becomes critical for clusters with low S/N such as this one. This cluster illustrates the relevance 
of point source subtraction when dealing with low number counts. In fact the point sources increase by about 60% the best-fit 
temperature of RDCSJ 1317-1-291 1 as one can see in Fig.|5] The low S/N implies also large errors on the luminosity and justifies the 
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disagreement between our results and ET04 shown in Fig. |4] The uncertainties on the luminosity are much larger than the difference 
expected in the luminosity when using different radii. 

4.5. ZWCL 1454.8+2233 

This cluster is very discrepant with respect to the behavior of other clusters in the L^ pz-T relationship given in Sect. |6] 
ZWCL 1454.8+2233 was identified as a relaxed cluster hosting a massive cooling flow by Allen et al] d 19961) using ASCA and 
RO SAT data. A 10 Ks Chandr a observation revealed the presence of two surface brightness edges on opposite sides of the X-ray 
peak which were discussed bv lMazzotta et alj (1200 lb under the hypothesis of a merging scenario. The 90 Ks Chandra observation 
analyzed by us confirms a very disturbed morphology. The surface brightness profile is inadequately described by a y6-model (the 
probability to accept the spatial fit is lower then 0.1%) both for the presence of the cooling core and for the presence of some surface 
brightness jumps. For this reason no values for Rc and /3 are indicated in Table|2l On the other hand for the spectral analysis a single- 
temperature model was accepted. The null-hypothesis probability is about 5%, a little less than the threshold indicated in Sect. 12.31 
We tried also to use a cooling flow spectral model mkcflow added to a mekal model but the improvement in the fit is minimal. The 
thermal complexity of this cluster, which could explain the peculiarity of its behavior in the Lhoi-T relation, convinced us to exclude 
it from the fit of the Lhoi-T relation. 

4.6. MS 1621. 5+2640 

The best-fit temperature that we obtained for MS 1621+2640 (kT - 7.5^q ' keV) is 10% greater than the temperature obtained by 
ET04 (kT = 6.8;^Q 5 keV) but in better agreement with the estimate by V102 (kT - 7.6^Q g keV). Also our luminosity extrapolated to 
R500 disagree with the estimate of ET04 by an amount on order of 18%. This discrepancy could be partly accounted for by the fact 
that we found a temperature of the gas higher with respect to ET04 temperature. 

4.7. CL 1641+4001 

For this cluster the effect of point sources on the cluster thermal emission is particularly important. In the region covered by the 
cluster we found six sources, two of which with a total flux > 50 x 10"'^ erg cm"^ s When the point sources are included in the 
fit the luminosity increases of ~ 50%. On the other hand these sources have a small (~ 3%) impact on the cluster temperature. 

4.8. MS 2053.7-0449 

MS 2053.7-0449 is assumed to be a cooling core cluster by V102. They measure a temperature of kT = 5.2;^^ !^ keV. Our best-fit 
temperature is ^T = 4.3^g4 keV within the Rspec region, and ^T = 5.1^} g keV when the larger radius Rgxt is used. Our lower 
temperature estimate could be explained by fact that the possible cooling region is not excluded in our analysis. ET04 did not 
exclude the cooling core and found a temperature of kT = 5.5^q 5 keV larger than VI02 and our temperature estimate. However, 
iMau^han et al. I (l2007l) found a value for the temperature of kT — 4.1^°-^ keV or kT - 4.2^^^ keV according if the central region is 
taken into account or not. From the comparison with the results of the above mentioned authors it seems that the cluster may not be 
a cooling core. 

5. Point source contribution to cluster thermal emission 

The effect of the point sources on the best-fit temperature is shown in Fig.|5] The plots show the estimated temperature for clusters 
plus point sources (Tduster+ps) versus the temperature for clusters without point sources (Tduster)- The left panel refers to temperatures 
estimated within Rspec while the right panel refers to Rext- Obviously the number of point sources within Rext is higher than within 
Rspec- The most discrepant point in Fig.|5]has been discussed in the note of RDCSJ 1317+291 1 (Sect.|4]i. 

The effect of the point sources on the temperature estimate is evident (especially within Rext)- In addition a mild dependence on 
the redshift is also present. The average of the ratio rx = Tdustei+ps/Tciusters gives a temperature excess of (8 + 3)% within Rspec and 
(13 + 4)% within Rext- SpHtting the sample at z = 0.7 we get, respectively, (16 + 7)% and (24 + 8)% for z > 0.7 and (3 + 1)% and 
(6+l)%forz<0.7. 

The left panel of Figure |6]represents the same plot of Fig. |5]for cluster luminosities: Lduster+ps is plotted versus Lduster- A different 
way of presenting these data is illustrated to the right of Fig. |6]where the ratio rL = Lduster+ps /Lduster ~ (1 + Lps/Lduster) is plotted 
versus Lduster- This is a way to examine the "contrast" of point source luminosity with respect to cluster luminosity. The plot shows 
a mild trend for the luminosity "contrast" to decrease with increasing cluster luminosity. Similarly to the temperature, the contrast 
in luminosity due to the point sources is more pronounced in distant clusters. The low-z clusters (open circles in Fig. |6] right) 
are generally below the high-z clusters (solid circles) in the same luminosity range. The average of rL gives an over-luminosity of 
(17 + 3)% with (22 + 3)% for high-z clusters and of (14 + 4)% for low-z clusters. These values are similar, but more significant, than 
those for rj computed within Rext- 

Given the probably different X-ray spectrum of clusters and point sources (possibly AGN) we have investigated the contribution 
of the point sources to the flux in the soft and hard energy bands separately, by computing the two ratios rs = Sduster+ps/Sduster for 
both 0.5-2.0 keV and 2.0-10.0 keV energy bands. The histograms of such ratios are shown, overplotted, in Fig.|2] where the solid 
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Fig. 5. Temperature of clusters plus point sources versus temperature of clusters without point sources. Left panel refers to Rspec and 
Right panel to Rgxt- Solid circles indicate high-z (z > 0.7) clusters and open circles indicate low-z (z < 0.7) clusters. In both panels 
the dashed line is equality between the two temperatures. 



(dashed) line represents the soft (hard) band. Clusters with redshift < 0.7 are shown to the left and those with redshift > 0.7 to the 
right. The figure indicates that the point source contribution is higher in the hard band, as one might expect if the point sources have 
a spectrum harder than clusters. This effect is more evident for the distant clusters (right panel). It has to be noted that in both energy 
bands the excess in flux is mainly due to some bright objects (projected onto or belonging to the cluster) rather than to many faint 
sources. 

The data suggest that there is a mild tendency for distant clusters (z > 0.7) to be more affected by point sources. If real, this 
effect would naively imply that the line of sight to farther clusters intercepts a higher number of bright point sources. However, 
this is contrary to what expected from the fact that distant clusters cover an average angular area which is ~ 50% that covered by 
nearby clusters, and that the flux is about half the flux of nearby clusters. If point sources are mostly unrelated to the clusters, these 
two effects would compensate and the point source contribution to the cluster overall budget should be independent of the cluster 
redshift. The observational result of seeing a larger number of point sources in distant clusters goes in the direction of having more 
sources belonging to high-z clusters. Indeed in a recent paper BR07 found an overdensity of point sources in clusters with respect 
to the field, and an indication in the hard band that the excess is mainly associated to high-z cluster. 



6. The Lio/-T relationship with and without point sources 

The effect of the point sources on L/,o/-T relation can be evaluated by direct comparison of the relations obtained by including 
and excluding the point sources. We first analyzed the Lio/-T using the gas temperature estimated within Rspec- Then, in order to 
consider the contribution on the temperature of all the sources within the extended region where the luminosity is computed, we 
used the gas temperature estimated within Rext- We express the Lio;-T relation as 

Lboi,44 - CTg (3) 
where Lboi,44 is the bolometric luminosity in units of 10"*"^ erg s"' and Te = T(keV)/6. 

The data and the fitted Lio/-T (using a.x^ method which takes into account the Lboi,44 and T uncertainties) are displayed in Fig. [8] 
for both the Rspec (left) and Rext (right), and the best-fit parameters are listed in Table [3] Solid and open symbols and solid and 
dashed lines represent quantities with and without point source inclusion. The most discrepant point in both plots is the cluster 
ZWCL 1454.8-H2233 which was excluded from the fit of the L^oz-T relation (see Notes on individual clusters in Sect. [Hi. The 
inclusion of the point sources, for both Rspec and Rext, does not have any significant effect on the slope and normalization of the 
Li)o/-T relation. This is because the correction to T and Lboi applied to each data point in the plot produces a shift in the T-Lboi 
plane almost parallel to the best-fit line. In Paper II (Branchesi et al. 2007) we analyze in more detail the Lm-T relation using data 
free of point sources. 
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Fig. 6. (Left panel) Luminosity of clusters plus point sources versus luminosity of clusters without point sources. The dashed line is 
equaUty between the two luminosities. (Right panel) Ratio between luminosity of the clusters plus point sources and luminosity of 
the clusters without point sources versus luminosity of the clusters alone. Solid circles refer to high-z clusters and open circles refer 
low-z clusters. 
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Fig. 7. Histograms of the ratio between the flux of the clusters plus point sources and the flux of clusters without point sources. The 
left panel is for clusters with redshift < 0.7 and the right panel for clusters with redshift> 0.7. The solid line indicates the soft energy 
band (0.5-2.0 keV) and the dashed line indicates the hard energy band (2.0-10.0 keV). 



7. Summary and Conclusions 

In this paper we have presented the details of the data analysis of a sample of 18 distant clusters (0.25 < z < 1.01) taken from the 
Chandra archive to derive the observational properties of the X-ray emitting gas. The same sample was used to study the point 
source counts in the inner region of distant clusters (BR07) and to study the evolution with redshift of the Lj,o;-T relation (Paper II, 
Branchesi et al. 2007). 
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Fig. 8. Cluster bolometric luminosity versus temperature relationship. The open circles indicate cluster plus point source emission 
and the solid circles indicate cluster emission alone. The dashed and solid lines are the best-fit of the Lh„i-T relation to the open and 
solid circles data respectively. The solid circle indicated as 1455 represents the cluster ZW CL 1454.8+2233 which was removed 
from the fit (see Notes on individual clusters in Sect.|4|l. 

Table 3. Lboi - T best-fit parameters 



Lbol - T(Rspcc) 


Cluster 


Cluster + Point Sources 


a 


+3.09-- 


+3 00 + " " 


logC 


+ 1.06 


+ 1 06 + "°" 


Xlijd-o.f. 


14.77/15 


17.54/15 


Lbol - T(Rexi) 


Cluster 


Cluster + Point Sources 


a 


+3 23 + 

-rj.^j - 0.38 


+3 29 + 


logC 


+ 1 04 + °'" 


+ 1 00 + ""^ 


Xljd-o.f. 


12.24/15 


14.25/15 



The very high angular resolution of Chandra that allows to isolate X-ray point sources embedded in the more extended X-ray 
emission from galaxy clusters, enabled us to estimate for the 18 clusters the effect of the point source non-thermal emission on the 
determination of the thermal emission due to the cluster itself. 

i) The point sources located within the cluster thermal emission region may affect considerably the estimates of X-ray observables 
like cluster temperature (by an amount up to 13%) and luminosity (by an amount of 17%). These percentages become larger if one 
considers clusters with z>0.7 where temperature and luminosity increase up to 24% and 22%, respectively (see Section |5]|. The 
results obtained suggest that, in order to estimate properly the observational parameters of the thermal emission, the point source 
contribution should be removed. 

ii) However, the inclusion (or exclusion) of point sources in the analysis of the hboi- T relation, indicates minor differences, within 
the uncertainties, of the relation (see Section |6l). This is due to the fact that for each cluster the correction to be applied to T and 
Lfeo/ for the presence of point sources produces a moderate shift in the Lio/-T plane almost parallel to the best-fit of the "correct" 
(excluding the point sources) L^oz-T relation. 
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Fig. A.l. {Left panel) The best-fit temperature obtained using the Cash statistics versus the best-fit temperature obtained with the x'' 
statistics. The dashed line is equality between the two temperatures. {Right panel) The luminosity obtained using the Cash statistics 
versus the best-fit luminosity obtained using the;^^^ statistics. The dashed line is equality between the two luminosities. 



Appendix A: Analysis of the best-fit parameters 

A.1. Comparison between Casli andx~ statistics 

In this Appendix we compare the best-fit temperatures and the bolometric luminosities obtained with the Cash and the statistics 
only for those spectrum files with the point sources excised. The agreement is remarkably good when the spectra are extracted 
within the region defined using Rspec (see Fig. lA.lb . The temperatures (luminosities) obtained with the Cash statistics are, on 
average, ~ 1% (~ 4%) higher then those obtained with the statistics. The histogram of the ratios between the temperatures 
(luminosities) obtained with thex^ and the Cash statistics is characterized by a standard deviation of 0.04 (0.03). 

The agreement remains good considering the more extended region defined when one uses the radius Rext- In this case the luminosi- 
ties with the Cash statistics are on average ~ 3% smaller than those obtained with^^, while the temperatures remain the same. The 
dispersion of the temperature histogram obtained within the region defined by Rgxt has a standard deviation of 0.1, higher than the 
dispersion obtained when using Rspec- This is presumably due to the fact that the spectral fit is better in the region that optimizes the 
signal to noise ratio. 



A.2. Best Parameters for Luminosity and Temperature 

The quality of the fits is better inside the region within Rspec, e.g the region that maximizes the signal to noise ratio. Thus the 
temperatures estimated within this region, Tr^^^ , are considered more representative of the actual average temperature of the gas. 
On the other hand the cluster luminosities are estimated within the extended region defined by the radius Rext in order to take into 
account the faint brightness tails at the cluster boundaries as well as the point sources in those regions. For each cluster the luminosity 
was calculated by fitting the counts accumulated within Rgxt with a thermal bremsstrahlung model. A best-fit temperature, indicated 
as Tr^|, is associated to each thermal model. Systematic differences between Tr^^^^ and Tr.^^ are within 3%. In order to understand 

how strongly the luminosity depends on the model and on its associated temperature, and in order to know which errors are involved 
when temperature and luminosity are calculated in two different regions (i.e. defined by Rspec and Rext, respectively) we computed 
the luminosity of each cluster in the Rext region but using the temperature obtained within Rspec- The results of this exercise are 
shown to the left of Fig. IA.2I The two methods give values consistent within the errors. The differences between the luminosities 
are all within 5%, that is within the statistical errors. 
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Fig. A. 2. Best-fit luminosity obtained within Rext with the temperature fixed to the value obtained within Rspec versus best-fit lumi- 
nosity obtained within Rext- The dashed Une is equality between the two best-fit luminosities. 



